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Abstract 

We present a unified theory of luminescence spectra for highly excited semiconduc- 
tors, which is applicable both to the electron-hole BCS state and to the exciton 
Bose-Einstein condensate. The crossover behavior between electron-hole BCS state 
and exciton Bose-Einstein condensate clearly manifests itself in the calculated lu- 
minescence spectra. The analysis is based on the Bethe-Salpeter equation combined 
with the generalized random-phase-approximation, which enables us to consider 
the multiple Coulomb scattering and the quantum fluctuation associated with the 
center-of-mass motion of electron-hole pairs. In the crossover regime, the calculated 
spectra are essentially different from results obtained by the BCS-like mean-field 
theory and the interacting Boson model. In particular, it is found that the broad 
spectrum, arising from the recombination of electron-hole BCS state, splits into 
the P- and P2-luminescence bands with decreasing the particle density. The depen- 
dence of these bands on the carrier density is in good agreement with experiments 
for highly excited semiconductors. 
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The issue of Bose-Einstein condensation (BEC) of excitons has received con- 
siderable attention because of the Bosonic nature of excitons [1]. Following the 
pioneering work by Keldysh et al. [2], the high-density electron- hole system 
with excitonic instability has been extensively studied both theoretically [3-5] 
and experimentally [6]. Recent developments in the experimental techniques 
make it possible to observe the remarkable optical phenomena suggesting the 
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spontaneous generation of the macroscopic quantum coherence. In particular, 
the anomalous exciton transport phenomena observed in CU2O [7] and Bil 3 
[8] are manifestations of the nondissipative rapid propagation of high-density 
excitons. However, current understanding of these experiments still remains 
controversial because of the complicated experimental situation due to the fi- 
nite life-time of excitons, phonon effects, spatial inhomogeneity of the system, 
and so forth. 

We should remark that the macroscopic quantum state of high-density ex- 
citons is qualitatively different from the BEC states of superfluid 4 He and 
atomic gases. The system exhibits a crossover from the exciton BEC to the 
electron-hole BCS state with increasing the excitation- light intensity. 

Under the weak photoexcitation, the system can be described by the inter- 
acting exciton model [3] because the mean distance between excitons is much 
larger than the exciton Bohr radius. The system undergoes the exciton BEC at 
sufficiently low temperatures, and the density of exciton with zero momentum 
is the order parameter of the macroscopic quantum state. When the photoex- 
citation is sufficiently strong, the bound electron-hole pairs can no longer be 
regarded as pure Bosons because the state-filling and the exchange effects, 
originated from the Fermionic nature of electrons and holes, take essential 
part in the problem. In the high-density limit, the electron-hole pairs behave 
like Cooper pairs, and the BCS-like energy gap at the Fermi level is the order 
parameter of the macroscopic quantum state [4,9,10] which is essentially the 
same as the excitonic phase [2]; it is noted that no Bosonic nature exists above 
T c , and the property of the system is totally different from the exciton BEC 
state. 

Similar crossover phenomena have been discussed in a variety of physical con- 
texts, e.g., superconductivity [11], nuclear matter [12] and superfluid 3 He [13]. 
Recently, much attention has been focused on the BCS-BEC crossover in con- 
nection with the unusual properties of the high-T c cuprate superconductors. 
In particular, it is observed that a pseudo-gap structure in the normal-state 
density of states of underdoped cuprates persists almost up to room tem- 
perature [14]. A further unusual property of the high-T c superconductor is 
its extreme short coherence length (of the order of several lattice constants) 
of Cooper pairs [15,16]; this fact is in contrast with the conventional super- 
conductors where the Cooper pairs are strongly overlapping in real space. 
When we investigate this kind of crossover phenomena, the optically gener- 
ated electron-hole system has the marked advantage over the superconducting 
materials, because the macroscopic quantum states can directly be controlled 
by changing the intensity or the frequency of excitation light. 

In this article, we first present a unified theory of luminescence spectra for 
highly excited electron-hole systems which is applicable both to the electron- 
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hole BCS state and to the exciton BEC state. The earlier theories based on 
the interacting-Boson model or on the BCS-like mean-field theory are not ap- 
propriate because the highly photoexcited system is often in the intermediate 
state between the electron-hole BCS state and the exciton BEC state. In the 
present theory, the multiple Coulomb interaction between electron-hole pairs is 
incorporated by numerically solving the Bethe-Salpeter (BS) equation [17] for 
electron-hole pair correlation function. This analysis is closely related to that 
in Ref. [18], where the absorption spectra for condensed exciton system are cal- 
culated by properly taking into account the ladder diagram for vertex function. 
In the present paper, the effect of the collective quantum fluctuation associated 
with the center-of-mass motion of electron-hole pairs is incorporated by the 
generalized random-phase approximation (RPA) [19]. The calculated spectra 
considerably deviate from results obtained by the BCS-like mean-field analy- 
sis, and are essentially different from results by the interacting Boson model. 
The crossover from the exciton BEC to the electron-hole BCS state is clearly 
found to manifest itself in the luminescence spectra. In particular, the P2- 
luminescence band, which arises from the exciton recombination accompanied 
by the excitation of another exciton from IS to 2S state, continuously changes 
to the Fermi-edge singularity with increasing the electron-hole density. In the 
low carrier density, we find that the intensity of the P 2 -luminescence band 
predominates comparing to that of the P-luminescence band, and this be- 
havior is in excellent agreement with experiments for various semiconductors 
[20,21]. The present analysis also predicts the weak luminescence band due to 
the recombination of electron-hole pairs which are generated by the collective 
phase fluctuation. These result are closely related to the recent noteworthy 
experiment of ZnO thin films, in which the room temperature ultraviolet laser 
emission is observed [22]. 

We consider a electron-hole system in a direct-gap semiconductor, which con- 
sists of the isotropic, nondegenerate parabolic conduction and valence bands 
with identical electron and hole masses. The repulsive interactions between 
electrons and between holes as well as the electron-hole attractive interaction 
are taken into account. The spin degrees of freedom is neglected to avoid the 
unnecessary complication due to the biexciton state. 

The spontaneous emission rate is expressed by the imaginary part of the 
electron- hole pair correlation function as I{ui) = — 2lmG(u — ji + 27), where 
7 and jj, are the exciton decay constant and the chemical potential of the 
electron-hole pairs, respectively. G(oo) is the Fourier transform of the electron- 
hole pair correlation function 

G{t) = -ie(t)Y,9k9p(ci(0)dl k (0)d_ p (t)c p (t)), (1) 

k,p 

where is the radiation- matter coupling constant; c p and d p are the anni- 



3 



hilation operators for electrons and holes, respectively. Here (• • •) indicates 
the expectation value with respect to the electron-hole BCS state. In order 
to consider the electron-hole pair correlation, let us introduce the Bogoliubov 
transformation given by c*, = u^k + Vk(3_ k , cL*, = UkP-k — VkOr k , where 
and P-k are the annihilation operators for new quasiparticles (Bogolons). The 
Bogoliubov parameters Uk and Vk are determined by solving the BCS-like gap 
equation for electron-hole systems. The BCS-like gap equation is reduced to 
the Wannier equation in the low density limit, which implies that the equation 
gives the correct result for the relative motion of electron-hole pairs both in 
the high- and low-density limits [11]. 

In the following analysis, it is convenient to introduce a two-component op- 
erator, ipk = (®k, P-k) an d the Pauli spin matrices (// = 1,2,3). The 
spontaneous emission rate is expressed as 



(w - yU) 2 + 7 2 4 



^2g k g v KkGk, P {u -fi + h)K p 

k,p 



(2) 



where K k = (u 2 k — v k , i). Here J coh is defined by 



^coh = Y,9k9pUkVkUpVp{(lpkT 3 1pk){lplT 3 'lp p )), (3) 
k,p 

and the (//, v) component of Gk,p(w) is defined by the Fourier transform of 



(Gk, P (t)) = -mml(o)r,Mo))(^p(t)ruMt))), (4) 

respectively. The first term on the right-hand side of Eq. (2) gives the the 
sharp spectrum at the chemical potential of the electron-hole pairs; this spec- 
tral component arises from the coherent spontaneous emission from the macro- 
scopic quantum state. The electron-hole pair recombination does not generate 
Bogolons so that the spectral width is merely determined by the lifetime of 
electrons in the conduction band. The second term expresses the incoher- 
ent spontaneous emission rate accompanied by the creation of Bogolons; this 
spectral component reflects the various interesting features arising from the 
many-body interaction in high-density electron-hole systems. 

Making use of the generalized RPA, we evaluate J co h in Eq. (2) as 



Icoh = \j2 9kU k Vk\ + 2^ 9k9pUkVkU p v p f ^4 vj_ ^ 

I fc J k, P _4 2m y z - h k - ti P ) 
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where V q (u) is the screened Coulomb potential. The single particle excitation 
energy (Bogolon energy) is given by E k = y^T+~A|, where ( k and A k are the 
renormalized energy of the electron-hole pair and the BCS gap, respectively. 
The first term on the right-hand side of Eq. (5) is obtained by the BCS-like 
mean field approximation, and the second term expresses the correction which 
comes from the collective quantum fluctuations similar to the Anderson mode 
in superconductivity [19]. 

The function Gk,p(w), which gives the intensity of the incoherent emission 
through Eq. (2), satisfies the following BS equation, 

g^(u)(ur + 2E kT2 ) - ]T v k _ k iu)g pM (u) (c^,r 2 + icj&n) 

k' 

= ((4(0)r^ fc (0))(^(0) T ^ p (0))), (6) 

where C k ° p = (u k u p + v k v p ) 2 and C k 3 p = (u k v p — v k u p ) 2 are the coherence 
factors and To is the 2x2 unit matrix. This BS equation is obtained by 
linearizing the equation-of- motion for the correlation function Gk,p(t)- The 
multiple Coulomb interaction between electrons and holes is expressed by the 
second term on the right-hand side of Eq. (6). This term makes a significant 
contribution to the recombination of the electron-hole pairs with nonzero total 
momentum, which is not taken into account in the BCS-like mean-field theory. 
The right-hand side of Eq. (6) is rewritten as 

((4(0)r^ fc (0))(^t( ) T ^ p (0))) 
= (r + r 2 )4, P - ir 3 E(-l) J+1 ((4(0)r^ p (0))(4(0)r^ p (0))). (7) 

Z 3=1 

The second term on the right-hand side represents the fluctuation effect due to 
the collective excitations similar to the Anderson mode in superconductivity. 
A generalized RPA analysis gives the following expression, 

3~ L — oo y 

(8) 

As a first step in our numerical calculations, the gap equation is iteratively 
solved to evaluate the renormalized energy of electron-hole pair ( k and the 
BCS gap function A fc . The incoherent luminescence spectrum is evaluated 
by numerically calculating the eigenvalues and eigenvectors of the stability 
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matrix [23] for Eq. (6). The quasi-static dielectric function given by the single- 
plasmon-pole approximation is employed [5,9]. 

The incoherent components of the luminescence spectra are shown in Figs. 1- 
3. As a reference, we also present the luminescence spectra calculated by the 
BCS-like mean-field theory and the generalized RPA analysis, in which the 
fluctuation effect associated with the center-of-mass motion of pairs is con- 
sidered but the multiple Coulomb interaction is neglected. We use the units 
where the exciton binding energy and the exciton Bohr radius being unity. 
The dimensionless mean interparticle distance r s = (47rn/3)~ 1//3 is introduced, 
where n is the electron-hole density. The exciton decay constant is chosen 
as 7 = 0.03, unless othewise stated. The arrow in the figures represents the 
level of the chemical potential of the electron-hole pair at which the coherent 
luminescence line appears. With decreasing the particle density, the position 
of coherent luminescence band approaches to uj — E g — — 1, where E g is the 
energy gap between the valence and conduction bands. This behavior ensures 
that the present analysis reproduces the result obtained by the interacting 
exciton model in the low density limit. 

Fig. 1 shows the luminescence spectrum in the high density case (r s = 1.0). 
The band renormalization effect due to the exchange interaction is clearly 
observed. The peaked structure near the Fermi level (oj — E g ~ 1) is the man- 
ifestation of the BCS-like gap formation, which is strongly pronounced by the 
multiple Coulomb interaction between electrons and holes; the electron-hole 
pairs with nonzero total momentum take part in the recombination processes. 
The luminescence band in u — E g > fi, which is missing in the mean-field 
analysis, corresponds to the pair annihilation process of the virtual Bogolons 
generated by the collective phase fluctuation. 

The luminescence spectra in the intermediate densities corresponding to r s = 
2.2 is shown in Fig. 2. It is noted that the incoherent luminescence spectrum 
below the chemical potential splits into two components: the broad component 
below uj — E g ~ —1.7 and the sharp one around oo — E g ~ —1.5. The intensity 
of the broad component is approximately proportional to the square of the 
particle density; its spectral position slightly shifts towards the shorter wave- 
length with decreasing the particle density. These behaviors suggest that the 
broad component corresponds to the P-luminescence band [20,21,24], which 
comes from the radiative recombination of an exciton assisted by the disso- 
ciation of another exciton. On the other hand, the sharp component around 
uj — E g ~ —1.5 corresponds to the P 2 -luminescence band [20,21], which arises 
from the radiative recombination of an exciton accompanied by the excitation 
of another exciton from IS to 2S state. In addition, a weak luminescence sub- 
band, called P3, is observed around u—E g ~ —1.6; this component comes from 
the exciton recombination assisted by the excitation of another exciton from 
IS to 3S state. The peaked structure above the chemical potential originates 
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from the pair annihilation of Bogolons, which are generated by the collective 
phase fluctuation from the electron-hole BCS state. It is noted that we do not 
deal with a simple two-exciton problem, but calculated P, P 2 , and P 3 bands 
in the framework of the many-body theory. We should also recall that these 
spectral components cannot be obtained by a BCS-like mean-field theory with 
the RPA correction. 

The luminescence spectrum in a low electron-hole density state (r s = 13) is 
shown in Fig. 3. It is found that the P-luminescence band almost disappears 
and the P 2 band predominates with decreasing the carrier density. This behav- 
ior arises because the transition of another exciton to the bound state (2S) is 
much stronger than that to the dissociated states in low carrier density state. 

We show in Fig, 4 the comparison between the present theory and the experi- 
ment for ZnO thin film at 77K [20]. Here the exciton decay constant is chosen 
as 7 = 0.1. The calculated spectra are in good agreements with this experi- 
mant in the following points: (1) the exciton luminescence line predominates 
in the low density state. (2) The intensity of the P 2 -luminescence line grows 
superlinearly with increasing carrier density, and it prevails in the interme- 
diate density states. (3) The P-luminescence line appears in higher densities 
and it merges with the P 2 -luminescence line as the carrier density increases. 

Finally, it should be noted that the dependence of the line shape and the 
spectral position on the electron-hole density cannot be predicted by the well- 
known analysis for two-exciton problem [25]. 

In conclusion, we first present the many-body theory for luminescence spectra 
which is applicable both to the exciton BEC state and to the electron-hole BCS 
state. Our analysis, based on the BS equation combined with the generalized 
RPA, properly describes the crossover behavior between the two essentially 
different macroscopic quantum states. 

This work is partially supported by a Grant-in-Aid for Scientific Research 
from the Ministry of Education, Science, Culture and Sports of Japan. 
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Fig. 1. The luminescence spectrum in the high-density state (r s = 1.0) calculated 
by (a) the BCS-like mean-field theory, (b) the generalized RPA, and (c) the present 
theory. 

Fig. 2. The luminescence spectrum in the intermediate-density state (r s = 2.2) 
calculated by (a) the BCS-like mean-field theory, (b) the generalized RPA, and (c) 
the present theory. 

Fig. 3. The luminescence spectrum in the low-density state (r s = 13) calculated by 

(a) the BCS-like mean-field theory, (b) the generalized RPA, and (c) the present 
theory. 

Fig. 4. The luminescence spectrum for various density states: (a) the present theory, 

(b) experiment for ZnO thin film at 77K (from Ref. [20]). Here Ex- line corresponds 
to the free exciton emission. 



9 



